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Abstract 

The role of external forces in systems exhibiting anomalous diffusion is 
discussed on the basis of the describing Langevin equations. Since there 
exist different possibilities to include the effect of an external field the 
concept of biasing and decoupled external fields is introduced. Comple- 
mentary to the recently established Langevin equations for anomalous 
diffusion in a time-dependent external force-field [Magdziarz et ai, Phys. 
Rev. Lett. 101, 210601 (2008)] the Langevin formulation of anomalous 
diffusion in a decoupled time-dependent force-field is derived. 

1 Introduction 

Over the last two decades it has become apparent that many complex systems 
exhibit a phenomenom which has been termed anomalous diffusion fl] El E] . On 
account of this, there has been an increasing interest in stochastic processes de- 
viating basically from standard diffusion processes characterized by a Gaussian 
behavior. Systems exhibiting anomalous diffusion differ from the linear time 
dependence of the second moment and rather show {x^) ^ t", where < a < 2. 
In this context processes with a > 1 dispersing faster than standard diffusion 
processes are called superdiffusive while a < 1 means that a system displays 
subdiffusive behavior. 

In the realm of anomalous diffusion, the classical diffusion equation has 
to be replaced by the so-called generalized diffusion equations [5]. The most 
prominent representatives of this class of equations are probably the fractional 
diffusion equations [3] , where the derivatives with respect to time or to space or 
both are replaced by non-integer order derivatives. A more fundamental account 
to anomalous diffusion is provided by a stochastic process called Continuous 
Time Random Walk (CTRW). This process generalizes the standard Random 
Walk and allows for random jump length and random waiting periods between 
the jumps [5]. It is well-known that the generalized diffusion equation can 
be derived from the governing equations of the CTRW. Another approach to 
anomalous diffusion has been put forward by Fogedby who proposed a coupled 
system of Langevin equations leading to the generalized diffusion equations . 
In a sense, this approach can be considered as a continuous realization of the 
CTRW. 



In the present paper we consider the effect of external forces onto processes 
exhibiting anomalous diffusion. Although the incorporation of external forces is 
straightforward in classical diffusion theory, leading to the well-known Fokker- 
Planck equations, this task appears to be rather involved when anomalous dif- 
fusion is considered. The arising difficulties are due to the long jumps and the 
long waiting times that can occur. Throughout this paper we distinguish be- 
tween biasing and decoupled external forces. This notation shall indicate that 
there are two different possibilities of the action of the force. When we speak 
of a biasing field, we mean that the external field acts as a bias only at the 
time of the actual jump. In contrast to this we speak of a decoupled field if the 
diffusing particle is affected permanently during the waiting time periods and 
hence the diffusion process is decoupled from the effect of the field. Note that 
this distinction is not necessary for classical diffusion processes. 

While the inclusion of external potentials is relatively well understood on the 
level of the generalized diffusion equations and the generalized Fokker-Planck 
equations respectively, there are still some open questions as long as the corre- 
sponding Langevin equations are considered. However, an exhaustive compre- 
hension of the Langevin equations is inevitable to investigate the properties of 
sample paths of such processes. 

The aim of this paper is to clarify the different possibilities of including an 
external force into the framework of anomalous diffusion, namely the difference 
between biasing and decoupled external forces, by considering the corresponding 
Langevin equations. It is organized as follows. First we state some fundamen- 
tals concerning the theory of anomalous diffusion and thereby shortly review 
Fogedby's continuous formulation of CTRWs and the concept of subordination. 
After shortly reviewing some wellknown and some very recent results on the 
Langevin formulations of the generalized Fokker-Planck equations for biasing ex- 
ternal fields we establish the Langevin equations for generalized Fokker-Planck 
equations for decoupled external potentials which have not been considered so 
far. We conclude with a discussion on the role of external forces in anomalous 
diffusion. 



2 CTRWs and Generalized Diffusion Equations 



A suited stochastic process to describe discrete sample realizations of many 
microscopic processes leading to anomalous diffusion is provided by the Contin- 
uous Time Random Walk (CTRW). This process is an extension of the standard 
random walk and allows for random waiting times between jumps of random 
length. In the decoupled case the CTRW is characterized by a waiting time 
distribution W{t) and a jump length distribution F{Ax). Depending on the 
properties of these distributions, the CTRW can lead to an anomalous behavior 
of the mean-squared-displacement [4j. The governing equation for the proba- 
bility distribution (pdf) of the position of the walker is the MontroU-Shlesinger 
master equation [7j 
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where the time kernel ^{t — t') is related to the waiting time distribution [5]. 
The master equation has a straightforward interpretation. It states that the 
density of particles at position x and time t is increased by particles that have 
been at x' at time t' and perform a jump from x' to x at time t. On the other 
hand the density is decreased by particles that have been at x and jump away 
at time t to some other position. The resulting process is non-Markovian for 
waiting time distributions with a power-law tail. 

Another account to describe the evolution of pdfs in the context of anomalous 
diffusion are the generalized diffusion or fractional diffusion equations. 

For a subdiffusive process the generalized diffusion equation can be cast into 
the form 

^^f(x, t) = dt'cj^it - t')D^f{x, t') . (2) 

For the special choice (/)(t) = (t)~^+", which corresponds to Mittag-Leffler type 
waiting time distributions W{t) ^ r~^~", Eq.([2]) yields the (time-) fractional 
diffusion equation 

|/(x,t)=P^"^/(.,t), (3) 

where T>1~°' is the Riemann-Liouville fractional derivative. It is well-known 
that generalized diffusion equations can be derived from the MontroU-Shlesinger 
equation, see e.g. [5]. In order to describe superdiffusive diffusion processes 
the so-called space- fractional diffusion equations have to be taken into account. 
These equations describe Markovian processes with power-law distributed jump- 
length and are often referred to as Levy flights. 



3 Fogedby's Approach and Subordination 

A continuous realization of the CTRW has been considered by Fogedby in [6]. 
His formulation is based on a system of coupled Langevin equations for the 
position X and time t 

x{s)^r{s), i{s)^Tj{s), (4) 

where T{s) and 77(5) are random noise sources which are assumed to be inde- 
pendent. ri(s) has in this context to be positive due to causality. The system 
Q can be interpreted as a standard Langevin equation in a internal time s that 
is subjected to a random time change. This random time change to the physical 
time t is described by the second equation. The combined process in physical 
time is then given according to x{t) — x[s{t)], where s{t) is the inverse process 
to t{s) defined as 

s{i) = inf{s : t{s) > i} (5) 

Closely related to this concept of Fogedby is the mathematical method of 
subordination. Using the (not to formal) notation of Fogedby one calls the 
process x{s) a parent process and s its operational time. The random time- 
transformation function t{s) has to be a non-decreasing right-continuous func- 
tion with an inverse function s{t). The resulting process in physical time t is 
then obtained by x{t) — x[s(t)] and is referred to as subordinated to the parent 
process. Consequently are the processes t{s) and s(t) named subordinator and 
inverse subordinator respectively. 
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In [3] it was shown that the Langevin equations ^ lead to a time-fractional 
diffusion equation if the 77(5) are governed by a generic one-sided a-stable dis- 
tribution. Generally the pdf of the subordinated process can be stated in the 
form 

f{x,t)^l dsp{s,t) fa{x,s) , (6) 
Jo 

where p{s, t) is the pdf of the inverse subordinator and fo{x, s) is the solution 
of the parent process [TUl Elj ■ 



4 Biasing External Force-fields 



Throughout this paper we will restrict to anomalous diffusion processes governed 
by waiting time distributions, i.e. we will consider equations of the form 
Hence we consider processes which are e.g. ruled by Levy-stable subordinators 
and time-fractional equations. The role of external potentials for Levy flights is 
discussed in [12]. 

Let us first clarify what we mean by biasing external forces. Therefore, 
consider the generic scenario of a subdiffusive CTRW governed by power-law 
distributed waiting times. A biasing external potential or force shall not affect 
the diffusing particle during the waiting periods but only provide it with a bias 
at the instance of a jump. In a sense one might say that the action of the force 
can be regarded as anomalous as well. 

If the considered force is time-independent it is well-known that anomalous 
diffusion in biasing fields can be described by the generalized Fokker-Planck 
equation 



d_ 

dt 



f{x,t)= f dt'<j>{t-t') 
Jo 
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where F{x) is the external force [31 [9]. The equivalent description based on 
Langevin equations is provided by the coupled system 



xis)^F{x)+T{s), i{s)^r,{s), 



(8) 



where T(s) is a Gaussian and 77(3) is a fully skewed a-stable Levy noise source 

If a time-dependent external force F{t) is considered it turns out that the 
situation is by far more involved. There exist different alternatives to include 
the force. One can for example consider the generalized Fokker-Planck equation 



^^f{x,t)= j\t'^{t-t') 



-^x^^'^^§x^ 



fix,t'). 



(9) 



However, such a generalized Fokker-Planck equation turns out to be physically 
meaningless. The correct equation has been found recently [HI HH \TE\ 

rt 



d_ 

dt 



f{x,t) = 



dx'^ 



dt'<i){t-t')f{x,t') 



(10) 



Recall that for the fractional time-kernel 0(t — t') — {t — t')^^' 
the fractional Fokker-Planck equation 



), Eq.dini) yields 



dx 



pi-"/(x,t). 



(11) 
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Figure 1: Time evolution of the probability density for a constant external 
biasing force F = v for Mittag-Leffler type waiting time distributions with 
a = 0.5 for the consecutive times t = 1 to 5. One can observe the persistence of 
the maximum at the origin indicating that there is no internal dynamic between 
the waiting times. The external biasing force acts only at the time of the 
displacements resulting in a asymmetry of the pdf. The plume stretches more 
and more into the direction of the force. 



Notice that the difficulty of time-dependent external forces stems from the fact 
that in this case the fractional derivative and the Fokker-Planck drift-term do 
not commute anymore. On the basis of the generalized Fokker-Planck equation, 
the difficulty is due the fact that it is not clear whether the external force 
has to depend on t' or t. For a detailed treatment of this issue, we refer the 
reader to the original papers. At this point, we want to confine ourselves to a 
simple plausibility argument to account for the correct operator ordering which 
naturally does not replace a rigorous derivation. 

As we have already mentioned, when time-dependent transition amplitudes 
F{x; x') are considered, the question arises whether this amplitude has to depend 
on t' or t which is equivalent to operator ordering problem in Eq. (fTO|) . To 
answer this question, let us consider the corresponding CTRW governed by 
the MontroU-Shlesinger equation ([T]). According to the interpretation of this 
equation the probability to be at the position x at time t is increased by the 
particles that jump at time t from some x' to x. Since this jump which occurs 
at time t is governed by the transition amplitude F{x\x') it is clear that the 
transition amplitude has to depend on the time of the jump, i.e. F{x]x' ,t). 
Performing the appropriate limit procedure, one obtains Eq. (jlOp as the correct 
FFPE for time-dependent Fokker-Planck operators. 

Consequently, the corresponding Langevin equation for a time-dependent 
forcing is not straight-forward to derive. In fact, it has even been stated in 
|14j that it is impossible to find a subordination description for time-dependent 
external fields. If the force is assumed to depend on the internal time, i.e. 

i(s)=F(s)H-r(s), Ks)^T^(s), (12) 

which corresponds to a completely subordinated force, the corresponding gener- 
alized Fokker-Planck equation would be Eq.® and thus Eq. fT^ lacks a physical 
interpretation. 
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The appropriate Langevin system has been found recently by Magdziarz and 
co-workers [TB]. They argued that a deterministic force should not be modified 
by the subordination procedure and depend on the physical time t and proposed 
the Langevin equations 

x{s)=F{t{s))+T{s), t(s)=77(s). (13) 

One recognizes that the force depends on the subordination process. Subor- 
dination of the process x{s) then yields for the force term i^[t(s[t])] = F{t) 
since t{s[i\) = t and hence the desired dependence on the physical time. It has 
been proven in [TB] that the Langevin equations (fT^ yields the same probability 
distributions as Eg.lfTU]) and hence that they describe the same process. 



5 Decoupled External Force-fields 

If a particle is assumed to be affected by an external potential throughout the 
whole waiting time period and the anomalous diffusion process is independent 
of this potential, we speak of a decoupled potential. 

It is instructive to consider a simple example where a particle which is ad- 
vected by a constant force during the waiting periods and performs jumps after 
the waiting periods. The pdf of such a process has been proven to be governed 

by 

fix, t)^ dt'cj^it - t')-^f{x - v{t - t'), <') , (14) 

which can be considered as a generalized advection-diffusion equation, where 
the advection is normal while the diffusion is anomalous [T7]. Observe the 
retardation of the pdf on the right-hand side which renders the equation non- 
local in space. A solution of this equation can be found after passing into a co- 
moving reference frame. The ansatz f{x,t) — F{^,t) with the shifted variable 
^ = X — vt yields a generalized diffusion equation for ^ 

|F(e, t) = 1^ dt'cj^it - t')^F{i, t') , (15) 

whose solution is given by (see Eq.Q) 

F{tt) = / dsp{s,t)Fo{^,s), (16) 

where Fq is the solution of the standard diffusion equation. 

In order to establish the corresponding set of Langevin equations, we have 
to be aware of the decoupled character of the advective field. That means, 
the advection has to be completely independent of the internal time s. Let us 
consider the following set of Langevin equations 

x(s) ^vri[s) +T{s), i{s)^r){s). (17) 

The solution of the subordinated process a;[s(t)] can be found by integration 

x{t) = x[sit)] = r^'^ V v{s')d[s'{t')] + B[s{t)] , (18) 
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Figure 2: Time evolution of the probability density for a decoupled external 
constant force with the same settings as in FigurelH i.e. Eg. ITU) with a fractional 
time-kernel. One can see that this evolution can be regarded as a force-free 
anomalous diffusion in a co-moving reference frame. Note that the force does 
not affect the diffusion process and hence can be considered as decoupled. 



where B[s(t)] means subordinated Brownian motion, that is the force-free pure 
subdiffusive part of the process [HI \W[ UHl [19] . The integral can be rewritten 
as 

rs{t) j-s{t) 

vri{s')d[s'{t')] ^ 

= f vdt' ^vt 
Jo 

(19) 

yielding for the subordinated process 

x{t) ^ vt + B[s{t)] . (20) 

Introducing the variable ^ = x — vt again, this equation can be written as 

m = B[sit)]. (21) 

Thus the variable ^ performs a force-free subdiffusive process and therefore 
yields the probability distributions given by Eq. (|16p . which proves that the 
Langevin equations ([T7|) actually corresponds to the generalized Fokker-Planck 
equation ([H|) . 

The case of time-dependent external field is only slightly more difficult. Con- 
sider a process, where the particle (of unit mass) performs during the waiting 
periods an overdamped motion according to the equation of motion 

x{t) = F{t) (22) 

where F{t) is some time-dependent force-field. The corresponding generalized 
Fokker-Planck equation reads [T7] 

f{x,t)^ j\t'<t>{t-t')^e~^>^'"^'''"^-f{xX)- (23) 
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The exponential function on the right-hand-side is the so-caUed Frobenius- 
Perron operator of the equation of motion for the deterministic part of x{t). 
This operator ensures the proper retardation of the probabihty distribution 
during the waiting period [20] , 

Since Eq. ((22|) describes an invertible conservative system Eg. ([23]) can be 
expressed as (see ^20] ) 

fix, t) = dt'cf>{t - t')^f{x - £ F{t")dt\t') . (24) 

Performing the ansatz /(x, t) ~ F{£,, t) with £, — x — F{t')dt' , the pdf of ^ is 
governed by the generahzed diffusion equation (|15p . 
The corresponding Langevin equation reads 

x{s) = F{s)T]{s) + ris), i(s)=r7(s). (25) 

Integration of this equation yields for the subordinated process 

.s(t) 

x{t)=x[s{t)] = / F{s')7^{s')d[s'{t')]+B[s{t)] 
Jo 

t 

F{t')dt' + B[s[t)] . (26) 

Evidently £, = x — J* F{t')dt' performs for this case a force-free subdiffusive 
process which proves that x{t) is a solution of Eq. fM)) . 

Note, however, at this point, that the inclusion of space-dependent forces is 
only straightforward as long as conservative dynamics is considered because only 
in that case the Frobenius-Perron operator can be expressed by a substitution 
operator like in Eq. (|24p . Even the simple case of a linearly damped motion 
between the random kicks, i.e. x — —^x leads to a generalized Fokker-Planck 
equation whose solution cannot be expressed in a closed form [17'. Hence the 
prove used here is not applicable anymore. Similarly, a closed form solution of 
the Langevin equation cannot be stated for this case. 

Comparing the Langevin equation for a biasing time-dependent force Eq. (|13p 
and the Langevin equation for the decoupled case, one realizes the difference 
between these equations. For the case of a biasing force, the force has to depend 
on the subordination process in the parent process. Then the force term yields 
the contribution F{t')ds{t') to the process. Observe that the force depends 
indeed on the physical time t but is integrated over the subordinated measure. 
In the decoupled case however, the force is integrated in physical time and thus 
is completely independent of the diffusion process. 

Of course it is possible to state the Langevin equation for a process where a 
biasing and decoupled force are acting independently. If Fb denotes the biasing 
and -Fb the decoupled force the corresponding Langevin equation reads 

x{s)^FD{s)7^{s)+FB{t{s))+T{s), Ks)^ij{s). (27) 

The time evolution of the pdf of such a process is displayed in Fig. [3] for a 
constant biasing force and a constant decoupled force with same amplitude but 
opposite sign. Note, however, that in many settings the two contributions are 
not independent and thus can display dependences. 
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Figure 3: Time evolution of the probability density for a decoupled external 
constant force combined with biasing constant force of same amplitude but 
opposite sign. Observe that the maximum of the distribution moves with the 
decoupled field while the distribution becomes more and more asymmetric due 
to the biasing field. 

6 Conclusions 

Concluding, in this paper we have discussed the effect of external forces on 
anomalous diffusion processes on basis of their corresponding Langevin equa- 
tions. We have introduced the concept of a biasing and a decoupled external 
field which has no classical analoge. Corresponding to the recently established 
Langevin formulation of biased diffusion in a time-dependent external field [16j , 
we have rigorously derived the Langevin equations for decoupled forces. To 
clarify the concept of biasing an decoupled external force in systems exhibiting 
anomalous diffusion we have presented the time evolution of probability den- 
sities for the different considered cases. We have shown that the established 
Langevin equation for decoupled force fields can be solved exactly for conserva- 
tive space-independent dynamics. 

The presented work has aimed at a clarification of the role of external forces 
in complex systems which are characterized by subdiffusion and long waiting 
times respectively. The approach based on the Langevin equation has provided 
thereby deep insight into the physical nature of the processes. 

Concluding we shall exemplify the concept by two simple applications each 
with a constant force. First consider the diffusion of tracer particles in an 
advective flow which has frequent obstacles such as e.g. sediments. In this case 
the external force, i.e. the advective flow, is decoupled from the diffusion process. 
Second, if the diffusion of groundwater through porous media is examined the 
gravity field provides a bias on the anomalous diffusion process. 
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